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ABSTRACT 



A surnmajry of the Welbull distribution and three repre 
sentative sampling procedures to determine point and con- 
fidence interval estimates of the parameters that occur 
in the functional form of the Welbull distribution are 
presented* Following this> a model, assuming the Welbull 
distribution, is proposed which could have possible ap- 
plications in analysing Polaris Missile System Trouble 
Failure Reports to determine point estimates of the 
reliability of the Polaris Missile components* 

I am indebted to Professor P. W* Zehna for his con- 
tinued patience, encouragement, and most capable guidance 
while acting as faculty advisor* I also wish to thank 
Professor W. M. Woods for his valuable assistance as 
second reader* 
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TABLE OP SYMBOLS 



(Listed In order of their use In the text) 
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Welbull distribution 
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the cumulative distribution function for a 
complex item 

the density function for a complex item 
the likelihood function 



V 



CHAPTER I 






INTRODUCTION 



This thesis Is concerned with sxxmmarlzlng the WelbuU 
distribution and providing several possible methods of 
analyzing Polaris Missile System Trouble Failure Reports 
to determine the reliability of Polaris Missile com- 
ponents when the Welbull distribution Is assumed to 
characterize the failure times of the components* 

The Welbull distribution, which was first proposed 
by Waloddl Welbull CiJ In 1951, is a member of a class of 
distributions which characterize wear out failure; the 
longer an Item has been used the greater the probability 
of failure* Wear out failure Is associated with gradual 
depletion of material, fatigue, or accumulated shocks 
and other factors* 

The key assiimptlon underlying many of the statistical 
proced\ires' In current use for evaluating the reliability 
of components or equipment Is that the failure times 
follow the exponential distribution* However, In many 
practical situations sufficient data Is not usually 
available to justify the assumption that the exponential 
failure law truly characterizes the time to failure 
distribution of the equipment or components being tested. 
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The Weibull distribution which includes the exponential 
distribution as a special case is a more general ass\imp- 
tion about the time to failure law of components or 
equipment* Kao has shown that the failure age 
distribution of electron tubes is better ohacterized 
by the Weibull distribution with shape parameter ap- 
proximately equal to 1.7 than by the exponential dis- 
tribution* Also when the exponential distribution is 
assumed to be the time to failure law* M* Zelen and 
M. C. Dannemlller [ 3 ] demonstrate that four commonly 
employed sampling procedures are not robust* with re- 
spect to Weibull alternatives* 

Chapter II summarizes the distribution theory 
associated with the Weibull distribution and defines 
terms to be used in the remainder of this thesis* 

Chapter III outlines sampling plans which can be em- 
ployed to obtain estimates of the parameters from ex- 
perimental' measurements* The various estimates to 

follow are derived under each of these plans. Chapter 

\ 

IV proposes three possible methods of using Polaris 
Missile Trouble Failure Reports for determining the 
reliability of Polaris system components* 

« Sampling procedures which are not very sensitive to 
departures from basic assumptions are called robust* 
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CHAPTER II 



SUMMARY OP THE WEIBULL DISTRIBUTION 



Lot X bo a random variable which may tako on any 
valuo in tho continuous scale from m to Infinity where 
m Is a non-nogatlve number. Tho three parameter Welbull 
cumulative distribution function (d.d.f.) Is defined as 



(1). P(x) 



a 

- 1-e for X t m, where 

0 X < m 



a s scale parameter, a > 0 
b - shape parameter, b > 0 
m X location parameter^, m % 0 
The probability density function Is defined as the 
first derivative of P(x) with respect to x, and is given by. 



(2). ' dP(x) - f(x) -- , 

dx . 



b(x-m) 



b-1 



-(x-m) 

a 



X m 



a 



V.0 



X < m 



Welbull fiJ demonst;pated that yield strength of 
Bofors steel, size distribution of fly ash, and fiber 



« Xq, m, and Xm, respectively were the notations used 
by Welbull 0L]] • 
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strength of Indian cotton are distributed in good agree- 
ment with the cumulative distribution function, equation 



function, equation (1), qualifies as a failure distribu- 
tion and has successfully employed the Weibull distribu- 
tion with location parameter, m, equal to zero to char- 
acterize the time to failure of electron tubes. For the 
purposes of this thesis and the sampling plans which 
follow in Chapter III, the location parameter, m, will 
be assvuned to be zero. This is a reasonable assumption 
since if an item is placed in operation at time zero, it 
is exposed to the risk of failure from the time it is 
first put into use. 

For completness of presentation, the moments and 
characteristic function of the three parameter Weibull 
distribution are listed below, but a detailed discussion 
will only be presented for the two parameter sub-family 
having msO, Given f(x) as the above probability den- 
sity function, the nth moment is defined by the formula 



(1), Kao has shown that the cumulative distribution 




-(x-m)'> 
b-1 a 



E(X*^) = 1 .3^^ b(x-m) 



e 



*<'**» a 

and is given by 



n 



k 




k 



A 

The mean and variance, denoted by P\ and <T respectively, 
are easily determined from the above relation to be 

1 

? I 1 

X s a j ' (b + 1) 

and 

+ 1) d + 1)J 

respectively. The function T( • ) is the gamma 
function which is defined for every v ?’0 by 

*“r— I C* 

~P(v) = ^ e“* dx 

Suppose now that T is a random variable which may 
take on any value In the continuous time scale frota zero 
to infinity, and the random variable T represents the 
total time to failure of an item. Let us further assume 
that the two parameter Weibull probability density func- 
tion, f(t), specifies the probability law of the random 
variable T where 

(3)« f(t) = (bt^**V®’) ®3cp (-t^/a), t ^ 0; a,b, >0 
— 0 elsewhere 

Accordingly, if an item is placed in operation at time 
zero, the probability that the total life of the item 
is less than or equal to t is expressed as 
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I"’ 





u 




I 



$mm 




P{T 5 t) = 




Therefor© the ounulatlve distribution function of the 
random variable T is given by 

P(t) = 1 - exp (-t'^/a) 

The reliability function^ denoted by R(t), of 
an item is also of particular Interest since it is 
desirable to know the probability that an item or 

I 

component does not fall in the time interval (0,t), 

Thus the reliability ftinction is the probability 
that the time to failure of an item is at least t 
time'' units and 1s expressed by 

R(t) 5 P(T > t) - (bt^“V«^)®3cp(-t^/a)dt, for t>0 

which upon integration becomes 
% 

( 5 ), R(t) 3 exp (-tVa) ; t 0 
Not© that the reliability function is merely the unity 
complement of the cumulative distribution f\motion» 
equation (If.)* 
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The nth moment of the random variable T is defined 
by the foirmla 

E(T“) s 5^t“(bt^“Va) o3cp (-tV®) dt 
and is given by 

n 

(6), E(Tn) = a J~* (f 1) 

The characteristic function of the random var- 
iable T is 

i\ 

c^(u) a E(ei'iT) s \ b t^*^ e“(tVa) dt 

which upon integration becomes 

C^(u) = ^ (lu) ^ “P 

The mean time to failure of the random variable T, 
denoted by X , is defined as the first moment about 
zero and Is easily determined from equation (6) to be 

1 

^ T **7 1 

(7). Xa:a I (S+1) 

In a similar manner, the variance denoted by 0"^of the 
random variable T is determined by the formula. 
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E(T^) - E(T)2, to bo 

(t -*• 1) -”P (F -f-1^ 

Another function of a time to failure probability 
distribution which provides a great deal of insight 
into the probability law is the hazard function or 
instantaneous failvire rate which will be denoted by 
Z(t)» The Instantaneous failure rate or hazard 
function is defined by ^ 

whore f(t) and R(t) are the density and reliability 
functions respectively. Substitution of f(t) and 
R(t) of equations (3) and (5) in the above express- 
ion, we have 

(8), Z(t) — bt^*V^ t>-0 

It is interesting to note that, for values of the shape 
parameter greater than one, the instantaneous failure 
rate is an increasing function of time, therefore indi- 
cating that wear out occurs during the life of the item. 
But in the special case where the shape parameter equals 
ohe and the Weibull distribution simplifies to the ex- 
ponential distribution, it is observed that the instan- 
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taneous failure rate Is a constant equal to 1/a, Inde- 
pendent of age. 

The effect of varying the shape and scale para- 
meters can most readily be illustrated by graphs of 
the probability density function because the prob- 
ability that the random variable T will lie in the 
interval (t, t + h) can be interpreted geometrically 
as the area under the probability density function in 
the interval from t to t-hh. Figure one Illustrates 
the probability density function for several values of 
the shape parameter, b, and with the scale parameter, 
a, equal to one. The effect of changing the value of 
the scale parameter is to merely squeeze or broaden 
the graph of the probability density function, l,e,, 
an Increase in the value of the scale parameter 
would broaden while a decrease in the value of the 
scale parameter would squeeze the graph of the prob- 
ability density fimction together. The first deriva- 
tive of f(t) with respect to t is 

f»(t) - b t^-2 ^ j".^ ^jb_^ 

which shows that f(t) has a maxlm\im at 
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if b is greater than one. But if the shape parameter 
has a value greater than zero but less than or equal 
to one, the probability density function is a mono- 
tonically decreasing function of t. Thus increas- 
ing values of the shape parameter squeezes the pro- 
bability density function together. This can be 
interpreted Intuitively to mean that the faster an 
item wears out, the greater the value of the shape 
parameter, b. Figure 2 Illustrates the reliability . 
function, R(t), for several values of the shape 
parameter and with the scale parameter, a, equal 
to one. Notice for fixed values of the scale para- 
meter, all the curves intersect at t equal to one, 
regardless of the value of the shape parameter. 

Thus far, consideration has only been given to a 
single item. It would be desirable to be able to con- 
sider a component or complex item consisting of M 
individuals or single items. Instead of assuming that 
the items are Independent, it will be assumed that the 
M items of the complex item are quasi- Independent in 
the following sense. If the life of the ith item is 
a random variable Tj^ with corresponding density func- 
tion Pj^(t) then let T be a random variable whose 
cumulative distribution function, G(t), is determined 
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R(t) 



' / 
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V 



Figure 2. Graph of reliability function for several values of shape parameter, b, and 

scale parameter, a, equal to one. 
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I 





by the equation 



l-G(t) = Pd^:^ t)P(T2:^ t) . . . P(Tj5>t) 



for t greater than or equal to zero 
the above equation for G(t) yields 



Solving 



T ^ 



t^i 

a< 



(9). 0(t) s 1-e , t^O 

The first derivative of 0(t) with respect to t de- 
noted by g(t) gives 



( 10 ). 



dG(t) ..g(t)s V '^l^ 

4r~^ 



bJ-l 




Equations (9) and (10) above will define the pro- 
bability density and cumulative distribution func- 
tion respectively of a complex ltem{^5j > l.e.* a 
component consisting of more than one item which are 

quasi-independent. The reliability function of a 

\ 

complex item is determined by the probability state- 
ment, P(T>t), to be 






xi 



(11). R(t)» \ g(t) dt=e , tP-0 

ahd the hazard function is determined from the re- 



lation Z(t).. g(t) to be 
l-G(t) 



) 
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(12). Z(t) = 

Thd formula for computing the nth moment of a com- 
plex item is then given by 

•OO 



E(T*') = ^ t^ (bj^t^^~Va^ ) exp 



M 

r 

'o Zl/ “ 

which can be integrated by parts to yield 



5 -® yVJ 

I nt”“^ exp (- ^ t'^i/a. )dt 

« ^ 

Again the mean and variance will be denoted by A and 
respectively and are determined by equation (13) to be 



(li^) 



r® M 



Jt C ^ w 

(15)* (T = j 2t exp (- 2! t”V®^i)<i^ * ^ 

Equations >(14) (1^) can be evaluated by nvunerical 

methods. 

If the Weibull family of probability distributions 
is assumed to characterize the time to failure of an 
item or component consisting of more than one item 
where the various members of the family differ only 
by the values of the parameters occuring in the func- 
tional form of the distribution^ it is clear that the 
reliability, hazard, density, and cumulative distri- 



but ion functions, along with the various moments are 
f\inotions of the parameters, a smd b. Consequently, 
useful probability statements about these items can 
only be made to the extent that estimates of these 
parameters may be given. At the present time, point 
and confidence inteinral estimation of parameters are 
most frequently employed. 

The maximum likelihood estimates of the para- 
meters, a and b, will be derived for each of the 
sampling procedures presented in chapters III and 
IV, although confidence interval estimates will only 
bo enumerated for the scale parameter, a, when the 
shape parameter, b, is asstimed to be known in one of 
the sampling plans. A thorough search of recent tech- 
nical publications indicates that confidence interval 
estimates, when both parameters are unknown, have not 
been derived and the derivation of these confidence 
interval estimates is beyond the scope of this thesis. 






CHAPTER III 
SAMPLING PLANS 

In thl? chapter, several methods of obtaining 
experimental data to determine the maximum likelihood 
estimates (M.L.E.) of the parameters that appear in the 
functional form of the two parameter Welbull probability 
distribution will be discussed. Experimental data can 
be obtained in a variety of ways, and each particular 
method of obtaining experimental measurements will be 
defined as a sampling plan. For each of the sampling 
plans presented, it is assumed that the random var- 
iable T characterizes the time to failure of an item 
and the density ftinction, f(t), defined by equation 
(3)» specifies the probability law of T. It is fur- 
ther assumed that it is possible to make repeated 
measurements, denoted by T 2 » ♦ • • , Tj^, of the 
random variable T, and the measurements, T]^, • . • ,Tjj, 
are independent identically distributed random var- 
iables, each with the same distribution as the ran- 
dom variable T# Accordingly, after selecting a ran- 
dom sample and performing an experiment, the observed 
values of the random variable T are designated by 

t^, ^ 2 * • • • 



) 
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is 






In the case of a complex Item composed of M 
individual items, a random sample of size n^ will 
be taken from the item of the 1th type, for 
1 = 1, 2, • • • ,M, Again the random variables 
Tij» for fixed 1 and j = 1, 2, , . . , n^,, are 
independent identically distributed random var- 
iables. Also after performing ah experiment, the 
observed values of the random variable will be 
denoted by t^2» • » '^l*^i* Each of the 

sampling plans presented are based upon random 
samples as defined in the above paragraph and are 
restatements of the sampling plans contained In 
reference , 

The method of maxlm\un likelihood is utilized 
to obtain estimates of the parameters, because this 
procedure generally tends to have desirable proper- 
ties not always shared by other methods. In brief, 
the method of maximum likelihood estimation consists 
of determining the values of the parameters that max- 
imize the likelihood function, denoted by L. The 
likelihood function, L, is a function of the para- 
meters in the joint density with fixed obseirvations. 



This maximizing process is generally accompolished by 
determining the values of the parameters that maximize 
the natural logarithm of L, since In L is an increas- 
ing function of L* The functional form of L will vary 
in general with each sarapl,ing plan, Lloyd said Llpow 



the maxlmiim likelihood method. 

The maximxim likelihood estimates of the para- 

A ^ 

meters, a and. b, will be denoted by a and b res- 
pectively in each of the sampling plans discussed 
in this paragraph. These sampling plans are out- 
lined as follows: 

SAMPLING PLAN I - Fixed sample size. 

Single item case - A random sample of size n is selec- 
ted and each member of the sample is tested until it 
fails. The total time to failure of each member of 
the sample is observed and recorded. If the obser- 
ved failure times of the random variable T are 
^1* ^2' • • • > likelihood fiinctlon 

is given by 



Upon substituting equation (3) for f(t^|^) in the 
above relation and taking the nat\iral logarithm 




and Mood [7] provide a detailed discussion of 



L= TT f(ti) 

is/ 
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of L, we have 



n n 

In L s n In b <- n In a + ^ In 

Therefore, If both the parameters, a and b, are un- 
known, the M.L.E. of those parameters are given by 

n A 

(16) , ss y“ t^Vn 

/s fy ^ ^ 

(17) . b ^ na/( ^ ^1^ In t4 - a Z In t.) 

JTt A,*» 

respectively, Lloyd and Lipow ^6^ outline a pro- 
cess of iteration which will solve equations (16) 

/ \ A A 

and (17) for a and b. In the event that the para- 
meter, b, is assvuned to be known, the M.L.E, of the 
parameter, a, under this sampling plan is given by 

(18) , ^ ^1^/n 

\ 

In this situation, it is also observed that a con- 
fidence Interval estimate for the , scale parsimeter 
can easily be determined using the statistic, 

(2 :E t>)/a 
* 

which has a Chi-square distribution with 2n degrees 
of freedom. Complete details of confidence inter- 
val estimates can be found in Mood E7l . 
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C omplex Item case - Assuming that the complex item is 
composed of M individual items, the procedure is to 
select a random sample of size from the item of 
the ith type, for Isl, 2, • * • ,M. Each member 
of the 1th random sample is tested until it fails, 
and the total time to failure is observed and re- 
corded for each member of the 1th sample. If the 
observed failure times of the random variable 
are tj 2 » • • • , ‘blUi, then the likelihood 

function for the item of the ith type is given by 



ni 






for i s 1, 2, • • • , M, Therefore, upon subs- 
tituting equation (3) for the above re- 

lation, we obtain 



Ljl = TT^^l exp (-tjLj^Vai) 

Accordingly, the M.L.E, of the scale and shape para- 
meters for the item of the ith type are determined, 
in the same manner as for a single item, to be 

(19), - .2 

and ^ 

r (20). 
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respectively, for 1 = 1, 2, . . • , M, Again, a 
process of iteration can be employed to solve for 



the majnimura likelihood estimates. In the situation 
where all of the shape parameters are assximed to be 
known, the M.L.E, of under this sampling plan is 
given by 



i*or 1 s 1, 2, • . . , M* 

sampling plan II - Fixed sample size with item 
tmmcation. 

Single item case - A random sample of size n is 
placed on te^it simultaneously, and the life test 
is terminated when r members of the sample fall. 
The n\amber of failures r is fixed in advance of 
the experiment, and r is an integer greater than 
zero but less than or equal to n. Using this 



sampling plan, the first r failures times of the 
random variable T are observed and recorded as 
the falltires occur. Then, if » 



Indicates that the likelihood function is given by 




i tp< oo are the first r failures, Halperin 
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where L la merely the density function of the first 
r order statistics. After substituting P(t) and f(t) 
of equation (3) and (i\.) in the above equation, the 
M,L,E, of the parameters, a and b, are detezmilned in 
the usual manner to be given by 




Again, the values of a and b which solve equations 
(22) and (23) can be determined by an iteration 
process. In the event that the shape parameter 
is assumed to be known, the M.L.E, of the para- 
meter, a, under this sampling plan is given by 

(24), » a ^ ( 2* ^1^ (n-r)tj,^)/r 

Complex item case - Assuming that the complex item 
is composed of H individual items, the procedure 
is to select a random sample of size n^ from the 
item of the ith type, for 1 > 1, 2, , • • , H, The 
n^ members of the ith random sample are placed on test 



1 
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simultaneously and the ith life test is terminated 
when members of the ith sample fall. As before, 
the number of failures r^ is fixed in advance of 
each experiment, and r^^ is an integer greater than 
zero but less than or equal to nj^* The first rj[ 
failures times of the random variable are 
observed and recorded as the failures ocotir. If 
the first ri failure times of the random variable 
Ti are given by 0 ^ t^i 2? t^2 - • • • 25: • 

then the likelihood function for the item of the 
1th type is 



/»JU 



ti - nil/(i»i - r^)l 7T 

V 






« 1 - **1 



for i’^1,2, Proceeding in exactly the 

same manner as with a single item, the M.L.E. of the 
' # 

parameters, a^^ and bj^, of the item of the ith type are 









A 

^1 






and 



aZ 



/s 

bl 



• ^26), ^>1 - 1 *♦* 



/V /l<Am 

(ni-ri) ^ 



IJ 
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As before, equations (2^) and ( 26 ) oan be solved by an 

A A 

iteration process for a^ and b^, for 1 =* 1, 2, , , , M, 
If all of the shape parameters are assumed to bo known, 
the M*L,E* of the parameter a^^ is 



SAMPLING PLAN III - Fixed sample size with time truncation. 
Single item case - A random sample of size n is selected 
and each member of the sample is tested for a specified 
time, say t^. The time t^, where t^ is greater than 
zero but less than infinity, is defined as the truncation 
time. The number of failures denoted by R that occur in 
the time intesrval, (0,tQ)^are observed and the R failure 
times of the random variable T recorded* Using this 
sampling procedure, R is of course a random variable 
which may take the value r, whore r s 0,1,2, * • * n. 

If the obsoiMTod failure times are t-^^jt^, • . • ,tp, 
then the likelihood function is 





for i = 1, 2 



M 




for t^ < t^ 



The likelihood function for this sampling plan is 
formally derived in Appendix A* Again, substituting 



P(t) and f(t) of equations (3) and (I 4 .) in the above 
relation, the M.L.E. of a and b are derived in the 
usual manner and determined to be given by 



respectively. As before, an iteration process can be 

A. /S 

employed to solve equations (28) and (29*) for a and b. 
Also, \mder this sampling plan the H.L.E. of the scale 
parameter is given by 



when the shape parameter is assumed to be known. 

Complex item case - If the complex item is assvuned to 

be composed of M individual items, the experimental 

procedure is to select a random sample of size nj^ from 

the item of the ith typo, for i = 1 , 2 , , , , ,M, The 

sample of the ith type is tested for a specified 

truncation time t4 ^and the number of failures denoted 

* ^o 

by that occur in the time interval^ (0,tj^^)j are ob- 
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served and the failure times of the random variable 
recorded. As before, the munber of failures, 
is a random variable which may take the value r^^, 

I 

where r^ st 1, 2, • « « , ^11 ^12* * * **^lr^ 

are the observed values of the r^^ failure times; then 
using the likelihood fxmctlon. 



the M,L.E, of aj^ and bj^ are detex*mlned to be given by 



t3l) 



( 32 ). 



+ (n,-r. )t^ In t4^ - a^ ^ In t 
1 1 ■‘■O •‘■O 4»l 



A A-w 

>■0 ^ 

\ 

for the item of the 1th type, where 1=1, 2, ,.,M. 
If all the shape parameters are assumed to be known, 
the M.L.E, of the parameter a^^ is 



Ij 



. ( 33 ). 





for 1 “ 1, 2, • • • # M, 



1 
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Thus, it is readily seen that considerable cal- 

A ^ 

culations are required to solve for a and b in each 
of the above sampling plans. But by utilizing a high 
speed electronic computer, the iteration process re- 
quired to solve for a and ^ is readily accompollshed. 
Kao ^91 outlines a graphical method of obtaining a 

A ^ 

first approximation for a and b. If these first ap- 

/V A 

proximations of a and b are provided to the computer 
as the Initial trial, the iteration method should 

A A 

quickly converge to a solution for a and b. In the 
next chapter, the above sampling plans, assuming the 
parameters, a suid b, are \mknown, will be utilized in 
the proposed models for analyzing Polaris Missile 
System Trouble Failure Reports, It will be further 
assumed that high speed electronic computers are 
available to aooompollsh the calculations required 

A A 

to solve for a and b. 
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CHAPTER IV 



APPLICATION TO POLARIS MISSILE SYSTEM 

In this chapter, three models are proposed which 
utilize Polaris Missile System (PMS) Trouble Failure 
Report (TPR) as experimental data for determining 
M.L.E. of parameters. By following the procedures 
proposed in the models below, it may bo possible 
to obtain a useful point estimate of quality indices, 
such as moan time to failure and reliability, of 
certain items composing the missile system. The 
models are applicable primarily to missile items 
for which a life history is maintained. An item 
will bo considered to have a life history if: 

(1) . It is possible to distinguish each item 
by some identification system from all other items 
of the same type, 

(2) , > A permanent record is maintained of the 
total operating time to failure, 

(3) * Assuming' that the item does not fail, it 

is possible to estimate accurately the total operating 
time accumulated by the item at some future date. 

In each of the models, the total operating time 
to failure of an item is assiuned to be a random variable . 
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^1 








T that is speolfied by tho two parameter Welbull prob- 
ability distribution with both parameters considered 
to be xinknown. Accordingly, If n Items of the same 
type are selected at random, this sample Is assiamed 
to consist of n Independent Identically distributed 
random variables which are characterized by the same 
probability law as T. It Is further assumed that 
none of the Items have accumulated operating time 
prior to their purchase by the United States Navy, 
Also, an Item, which has failed and been repaired. 

Is considered to be the same as a newly purchased 
Item with no accumulated operating time. In each 
of the models, the operating time of the Items Is 
accumulated by check outs or other reasons for op- 
erating the Item, and the TPR which provides the 
total operating time to failure Is used as the ob- 
served value of the random variable. Naturally, it 
is assxamed' that the time to failure reported in a 
TPR is accurate. 

The first model proposed Is Sampling Plan .1 
which is discussed In Chapter III. Consequently, 
for a sample of size n, the M.L.E, of a and b are 
given by equations (16) and (17)» respectively, 
and a point estimate for the reliability and mean 
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time to failure can be obtained by substituting the 

A 

calculated values of a and b into equations (5) and 
( 7 )» respectively* In the case of a complex item, 
the procedure is the same, except that the values 
are obtained for the M.L.E, of a^^ and bj^ for 
1 a 1, 2, • • * , M, with equations (19) and (20), 
Again, a point estimate of the' reliability and mean 
time to failure can be obtained by substituting the 
values of a^^ and bj^ into equations (11) and (II4.), 
respectively. Notice that equation (ll^.) must be ; 
evaluated by numerical methods to obtain a point 

f 

estimate of ^ • The limitation of this model is 
the waiting time required to observe the nth failure. 
Prior to outlining the next two models which 
are applications of Sampling Plans II and III, it is . 
necessary to present additional background concemlng 
the PMS, Consider an item of a given type, say X, 
which continues to accumulate operating time by 
check outs \intil it falls* Because the PMS is an 
expanding program, items of type X are constantly 
being added to the system and also acciunulating 
operating time. Consequently, if n items of t3rpe 
X are selected as a sample, each of the items could 
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have accumulated varying amounts of operating time 
depending upon when the item was purchased by the 
Navy, The effect of the above situation is to com- 
plicate Item and time truncation sampling plans as 
discussed In Chapter III. But Sampling Plans II and 
III can still be used to obtain M.L.E. of the para- 
meters ,lf a life history Is kept for Items of type 
X» as follows. 

The second method for obtaining point estimates 
of quality indices employs Sampling Plan II, The 
procedure is to select P samples of size nj, where 
j s 1, 2, , , . ,P, from items of type X, The items 
In the Jth sample are selected so that each of the 
members of this sample are expected to have accumu- 
lated approximately the same amount of operating 
time at a future date. For each of the P samples, 
fix in advance the number of observed failures, de- 
noted by rj, which will terminate the jth experiment. 
If equation (21) Is used as the likelihood function 
for the jth sample of size n^ and if the P samples 
are assximed to be Independent, the likelihood 
function, denoted by ^ , of the P samples is given by 
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I 





Upon substituting P(t) and f(t) of equations (3) and 
(14.) in the above relation, the M.L.E. of the para- 
meters, a and b, are determined in the usual manner 
to be given by 



Therefore, a point estimate of a quality index can be 
obtained by substituting equations (31+) and (35) into 
the equation which defines the quality index desired. 

In the case of a complex item, the method is completely 
analogous except there are M items of different types. 
Again, life histories must be maintained for each type 
of item. ' 

The third model uses Sampling Plan III. As in • 
the second model, P samples of size n^, where 
j -1, 2, • • • ,P, are selected from items of type 
X so that each member of the Jth sarple is expected 
to have acctimulated the same amount of operating 
time at a future date. A truncation time, denoted 




7 
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by is selected for each of the j samples and the 

experiment is considered to be terminated after the 

truncation time has elasped for each of the -j samples 

of size n^, where j -1, 2, • • * ,P. As in Sampling 

Plan III, the number of failures, denoted by Rj, which 

occur during the truncation time, t < , in the jth 

•'o 

sample is a random variable. Assuming that each of the 
P samples are independent, the likelihood function, de- 
noted by , is 

• 4 r-^L, 



Where is the likelihood function of Sampling 
Plan III, Proceeding as in Sampling Plan III, 
the M,L,E, of a' and b are given by 




Gihd 
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In the case of a complex item the procedure is com- 
pletely analogous. Therefore, point estimates of 
quality indices can be determined in both the single 
and complex item ease by substituting the equations 
for the M.L.E, of a and b into the equation for the 
desired quality index. 

In the last two models, the selection of the 
members of the jth sample should occur naturally 
since a group of items of a given type is usually 
issued to a submarine at the same time. These 
items should acc\amulate operating time at the seime 
rate through check outs \antil failure. The time 
when each of the items is issued can be detemined 
by examining the life history of the items; con- 
sequently, the members of the jth sample can be 
determined for either model two or three. 
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APPENDIX A 



In Sampling Plan III, n items are tested vmtil a 
specified time, t^, has elasped. The number of failures 
denoted by R that occur in the time interval (0,tQ) are 
observed and the R failure times recorded. Using this 
sampling procedure, R is a random variable which may 
take the value, r, where r = 0, 1, . , , , n. As 
before, equations (3) and ik.) define the probability 
density and cumulative distribution fvuictions res- 
pectively of the random variable T. The discussion 
which follows is a formal derivative of the likelihood 
ftinctlon for Sampling Plan III. 

Maximum likelihood estimators are functions of 
the data obtained from an experiment. In this case, 
the experiment is to specify a truncation time t^ 
and observe the number of failures, R, and the 
failure times, Tj^, which occxir in the time in- 
terval (0,tQ); consequently, the joint density 
function of the number of failures, R, and the 
times to failure, given that the failures 
occ\ir in the time lntex*val would bo ap- 

propriate for determining the likelihood function. 
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Let A be defined as the event, 

^2^ ^2» • • • t ■fcjj) where the Tj^ 
are not ordered and t^< t^ for i=l, 2, . . ,,R. 

Let B be defined as the event, (Rsr), where 

3^- 0, 1, . , , , n. 

Let C be defined as the event that the failures 
which occxxr do so before time t^. 

Let the event, AB, be expressed as t^, . . . ,Tj, < t^) 

Then the problem is to obtain the distribution ftinctlon 
which represents the probability statement P(AbIc). 

Using conditional probability relations, it is ob- 
served that 

P(ABlO- P(ABC) . P(BC) _ P(AlBC>P(B|C) 

P(BC) p(c) " 

where 

?U\BG)^ P(Ti<ti, . . . ,TR<tR|BnO 

* * * »^r*^^r> Jr-fi^^o* » » » >T n>to}nB/T£j 

P(Bnc) 
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and 



p(b|c) = jJ-P(to^ 

therefore 

F(ti,t2,...,tr,R=r)= |J-P(to)J 

for ti^^^ 

Accordingly, the joint density function is obtained 
by taking the partial derivatives of the above dis- 
tribution function with respect to tj: for isl,*..,r 
and is given formally by 



f (t^, • • • , tjt,R s r) 



for t t^ 
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